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The dynamics of a coupled two-component nonequilibrium system is examined by means of 
continuum field theory representing the corresponding master equation. Particles of species A may 
perform hopping processes only when particles of different type B are present in their environment. 
Species B is subject to diffusion-limited reactions. If the density of B particles attains a finite 
asymptotic value (active state), the A species displays normal diffusion. On the other hand, if the 
B density decays algebraically oc t'" at long times (inactive state), the efi^ective attractive A-B 
interaction is weakened. The combination of B decay and activated A hopping processes gives rise 
to anomalous diffusion, with mean-square displacement {xA{t)^) oc for a < 1. Such algebraic 
subdiffusive behavior ensues for n-th order B annihilation reactions {nB 0) with n > 3, and 
n — 2 for d < 2. The mean-square displacement of the A particles grows only logarithmically 
with time in the case of B pair annihilation (n — 2) and d > 2 dimensions. For radioactive B 
decay (n = 1), the A particles remain localized. If the A particles may hop spontaneously as well, 
or if additional random forces are present, the A-B coupling becomes irrelevant, and conventional 
diffusion is recovered in the long-time limit. 
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I. INTRODUCTION 

There has been considerable effort to elucidate the 
properties and conditions of anomalous diffusive behav- 
ior. A simple physical realization is given by diffusion 
on a fractal lattice , where due to the increasing num- 
ber of paths within the lattice, the time for a diffusion 
process will be prolongated. Also diffusion in random 
media with quenched disorder may be anomalous. De- 
pending on the distribution of barrier heights (or depths 
of traps) one may observe normal diffusive or subdiffu- 
sive behavior, respectively, if the available number of dif- 
fusive paths is reduced by the presence of obstacles [^,^ . 
Here we discuss a quite different situation in which dif- 
fusion is activated by the presence of particles or excita- 
tions which also propagate diffusively, but in the course 
of time decay. As a result the activated diffusion is ren- 
dered anomalous because the number of available paths 
decreases with time. However, the resulting structure 
of diffusive paths is not static, but evolves temporally. 
One may call this phenomenon dynamical fractality or 
dynamical disorder, depending on how the spatial distri- 
bution of excitations evolves in time. 

We model this scenario by starting from a two- 
component system consisting of distinct particle species 
A and B, with local time-dependent densities pA{x,t) 
and pBix,t). An A particle is allowed to perform hopping 
processes between adjacent neighboring sites on a lattice, 
provided there are one or more B particles present in its 
vicinity. To be more specific, an A particle hops from 
a site j to a neighboring point i subject to the condi- 
tion that this site i is already occupied with a particle 
of species B, and with a rate proportional to the local B 



particle density. Obviously, such an effective attractive 
interaction strongly influences the diffusive mobility of 
the A species: Their mean-square displacement {xA{t)'^) 
will depend on the time evolution of the local B density 
Pb{x, t). 

A non-trivial temporal behavior for ps {x, t) will result 
if we submit the B species locally to diffusion-limited 
reactions such as n-th order annihilation nB ~^ (at 
the same or adjacent lattice points) or combined anni- 
hilation (n > 2) and spontaneous offspring production 
B ^ {rn + 1)B (the B particles then perform branch- 
ing and annihilating random walks, BARW). Once the 
time-dependence of psit) has been determined, we shall 
see that the A kinetics is in the long-time limit to good 
approximation described on the basis of the associated 
mean-field rate equation. When the B species is in an 
active state, i.e., pB{x,t — *■ 00) = > 0, with a basi- 
cally homogeneous distribution in space, the A particles 
will display normal diffusive behavior, with a diffusion 
constant Da oc p^. In such a situation one has dynami- 
cal disorder, but there is always a finite fraction of sites 
available for hopping. However, an inactive phase, or the 
BARW critical point, are described either by an expo- 
nential decay pB{x,t 00) oc e~^*, in which case the A 
particles remain localized, or by a power-law decrease 
PB{x,t — > 00) oc t~°' with a characteristic exponent 
a > 0. The diminishing density of B particles reduces 
the induced mobility of the A species, and these compet- 
ing effects lead to subdiffusive behavior {xA{tY) oc t^~" 
for a < 1. In the borderline case a = 1 one has merely 
logarithmic growth {xAit)"^) oc Ini. 

Intuitively some of this behavior can be easily under- 
stood in the limit of vanishing diffusion of the B-particles. 
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We allow for multiple occupation, i.e., the site occupa- 
tion number can be any integer between zero and infinity. 
When the occupation number of a i? particles at a cer- 
tain lattice point is nonzero, that site is available for the 
A species. In an inactive state the number of B parti- 
cles will be permanently reduced by reactions leading to 
a decreasing density of available sites for the A species. 
This procedure can be viewed as an effective "thinning 
out" of lattice sites that leads to subdiffusive behavior 
for the A particles reminiscent of (but distinct to) the 
mathematically considerably more complex phenomenon 
of diffusion on a fractal lattice. 



II. THE MODEL 

Here we present a more precise definition of our model 
in terms of a master equation which we formulate in the 
standard Fock-space formulation £ ] , sometimes called 
"quantum Hamiltonian formalism" |7,||, particularly for 
particles with hard-core repulsion. Physically, this inter- 
action is usually insignificant unless exclusion between 
different particle species or external driving forces need 
to be taken into account. This is intuitively clear for an- 
nihilating processes where the particle density tends to 
zero at late times [p[-p^, but remains true also in the 
absence of particle reactions even in one dimension. 
In some models, however, e.g., the annihilation-fission 
model 0] or pair contact process with diffusion |p^ , |l6| , 
site occupation number constraints do play a crucial role. 
Hence, usually the specific choice of model is prescribed 
by the mathematical treatment used to analyze it. In the 
present context where we shall employ mean-field tech- 
niques and renormalization group arguments it is more 
advantageous to consider particles without site exclusion 
(for recent reviews, see e.g. Refs. 0,0). 

We consider a system consisting of two different types 
of particles denoted as A and B. The time evolution can 
be represented through an evolution operator L 
The corresponding annihilation and creation operators 
are written as (bi) and aj (bl), where the index indi- 
cates a lattice point in d space dimensions. For example, 
the normal hopping process of species B from a site j 
to its neighbor i is described by the evolution operator 
D {blbj — b^j), and for the entire lattice therefore 



Lr = D 



where D is the hopping rate or diffusion constant. 

An analogous expression would describe free diffusion 
of the A particles. Here, however, we examine the situa- 
tion that such a process is only allowed if there is at least 
one B particle present at site i. If no representative of 
species B is available at that site, an A particle cannot 
move there. The time evolution operator for that process 
is proportional to (a^a^ — ajaj)6|^5i. The corresponding 



hopping process will occur provided an A particle is in 
fact present at site j and at least one B particle occu- 
pies site i. Moreover, its rate is actually proportional to 
the number of B particles present at site i. For the full 
system we obtain 



D 



4b^ 



ai fo] bj 



(2.2) 



Here, D denotes the induced transition rate for the dy- 
namical process of species A. 

In contrast to species A, the B particles are subject to 
local reactions. A decreasing number of B particles will 
lead to a slowing down for the motion of A's through the 
lattice. For n-th order annihilation reactions nB 0, 
the nonequilibrium evolution operator reads ||l^ 



Li. = A„5:(i-C)^r. 



(2.3) 



Obviously, this operator describes the annihiliation of 
n particles of type B at a, lattice site i provided such 
particles are available; A„ denotes the corresponding 
rate. Similarly, spontaneous branching processes B 
(to + 1)B with rate a„i are described by p9| 



Lp = (Tr 



i:('>r 
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b]b. 



(2.4) 



Together with Eq. (p.l|), Lr and Lp represent the time 
evolution operator for branching annihilating random 
walks (BARW). 

The complete dynamics is determined by i = La + 
Lb + Lfi(+Lp), and may be encoded into a time- 
dependent "state vector" 



|^^(t))^^P(n,t)|n) 



(2.5) 



Here P{n, t) is the evolving probability distribution for 
the unrestricted site occupation numbers n — {rii} for 
both A and B particles, and \n) a basic vector containing 
all possible entries rii ~ 0, 1, 2, ...oo, i.e., the eigenvalues 
of the second-quantized bosonic particle number opera- 
tors a|ai and bjbi, respectively. The state |0) represents 
the vacuum with no particles present, |0) = = bi |0). 
The state vector obeys the equation of motion 



dt\F{t))=L\F{t)) 



(2.6) 



(2.1) or formally \F{t)) 



The nonequilibrium operator L corresponds to, and is 
obtained from the evolution operator L' of the classical 
master equation that can generally be written as 



dt P{n, t) = L' P{n, t) , 



(2.7) 



and the matrix elements of L and L' are uniquely related 
to each other. The time-dependent average of an arbi- 
trary physical quantity G{n) with the probability distri- 
bution P{n, t) can be cast into a "matrix element" form 
for the corresponding second-quantized operator G(t) 
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(G(i)) = E ^("' = (*l 1^^ W) , (2.8) 

with the projection state = (0| exp^^(ai + 6;). Us- 
ing the relation {'i'lL ~ 0, the evolution equation for an 
arbitrary operator G becomes 



d,{G)^{^\[G,L]\F{t)) 



(2.9) 



All the dynamical equations governing the classical prob- 
lem are thus determined by the commutation rules of the 
underlying operators and the structure of the evolution 
operator L. In our case the dynamics of the model is 
given by induced hopping processes for the A particles 
and diffusion-limited reactions for the B species, which 
we shall assume to be distributed randomly at the initial 
time t = 0. 

As a final step, we emplo y co herent basis states to rep- 
resent the matrix element ( p.8| ) by means of a path inte- 
gral PJl7||, and take the continuum limit. Absorbing fac- 
tors containing the lattice constant into the diffusion and 
reaction rates, we may compute averages with a dynam- 
ical weight exp(— y^[a, a, b, 6]) consisting of contributions 
to the bosonic field action A describing the ordinary B 
diffusion 



Asib, b]^ J d'^x jdtb {dtb - D \/^b) 
the pure n-th order annihilation reactions 



(2.10) 



Anib, b] = -A„ / d'^x dt{l~ 6") 6" , (2.11) 



or offspring production processes, 

Ap[b, b] = cfm I d'^x [dt(l- S™ ] bb , 



(2.12) 



respectively. Finally the A diffusion as induced by the 
coupling to the B species. 



^yi[a, a, b,b] — / d 



dt( 



dta~ D{V^a)bb 

+DaV^{bb)\ . (2.13) 



Notice that b{x, t) b{x, t) represents the local density 
Pb (a;, t) (when appropriate ensemble averages are taken) ; 
the A diffusion is thus mediated by the presence of B par- 
ticles. We remark that apart from the continuum limit, 
the mapping of the master equation onto the above field 
theory is exact and involves no further approximations. 
(We have omitted the boundary contributions stemming 
from the initial conditions and the projection state here.) 



III. GENERAL CONSIDERATIONS 

Clearly, the dynamic process for the A particles as de- 
fined above is induced by the coupling to the reactive 



B species only. When there are no B particles present, 
PB{x,t) = 0, the A dynamics obviously ceases. Indeed, 
it turns out that there appears no noise in the dynamic 
equation governing the A kinetics, which would formally 
appear as a contribution oc a a (or higher powers of d) 
in the dynamic functional. In fact, any stochasticity 
emerges as a result of spatio-temporal fluctuations for 
the B species (essentially reaction noise here). In order 
to further elucidate this point, we may derive effective 
Langevin-type equations for the local densities pA and 
Pb- To this end, we need to perform the shifts d = 1 -I- a, 
b — 1 + fe, which take care of the annihilation operators 
appearing in the projection state (^'|, see Ref. ||l^. To 
be specific, let us consider the case of B pair annihilation 
reactions. Omitting temporal boundary terms describing 
the initial configuration, the new action becomes 

y^[a,a,6,6] = jd'^xjdt a(dta - D' V^a ~ D (V^a) b 
+Da{V^b)) - Da{V^a)bb + DaaV^ibb) 



-b {dtb - DV^b + 2A b^) + X b^b^ 



(3.1) 



Here, we have allowed for additional ordinary A diffusion 
processes with rate D' . This dynamic action is equivalent 
to the following set of coupled Langevin equations, 

dta ^ D' V^a + D (V^a) b~Da (V^fe) + C , (3.2) 
dtb^ D\7^b-2Xb'^ + r] , (3.3) 

where the fiuctuating forces with zero mean are charac- 
terized by the noise correlations 

(C(af,t)C(^',i')) =0^, 

{({x, t) Ti{x,t')) = b [V'^a{x, t)] b{x, t) 5{x - x) 5{t - t') 
-b a{x, t) [6(x, t) 5{x - x) 5{t ~ t')] , 
{f]{x, t) 7j{x', t')) = -2 A b{x, tf 6{x - x') 5{t ~ t') . (3.4) 

Taking averages, we may then identify paU) — {o-{x,t)) 
and ps(t) = {b{x,t)), as Eqs. (3.2) and ( p^ ) obviously 
generalize the mean-field rate equations for the local par- 
ticle densities. The reaction noise for the B species 
displays the characteristic negative correlations ("imag- 
inary noise"), which reflect the particle anticorrelations 
induced by the annihilation reaction |18 2C| , p^ , [l^ . When 
there are no B particles left [b{x, t) — Q], the fluctuations 
cease, characteristic of an absorbing inactive state. As 
anticipated, no noise contributions exist for the pure A 
dynamics, but there appear A-B noise cross-correlations. 
(Notice that pure diffusion noise does not appear explic- 
itly here.) 

Next, let us study what happens when the A particles 
are subject to an additional random force that leads to 
ordinary diffusion, i.e., the term oc D' in the action ( [3.l| ). 
Obviously, one should expect that the induced diffusion 
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cx D is suppressed in this situation, and in the long-time 
hmit standard diffusion prevails. This becomes indeed 
clear through simple power counting, introducing a mo- 
mentum scale K, i.e., [x\ — k~^, and measuring time 
scales as [t] — k^'^, as appropriate for diffusive dynamics. 
Then [D] = [D'] = become dimensionless, and we infer 
the field scaling dimensions [a] = [6] = [a] = \b] = k° and 



as to be expected for d-dimensional par- 



ticle densities. The remaining couplings (reaction rates) 
acquire the scaling dimensions 



[A„] = K 



^ 2-{n-l)d 



[D] = 



(3.5) 



A positive scaling dimension means that the corre- 
sponding parameter is relevant in the renormalization- 
group (RG) sense. E.g., the branching rate Um carries 
the dimensions of a "mass" term, and indeed represents 
the decisive control parameter for BARW: In mean-field 
theory, the critical point must be at (T,„ = 0, and is there- 
fore described by the pure annihilation model, while for 
any positive am there will be only an active phase char- 
acterized by exponential correlations. The annihilation 
rate is relevant for d < 2/ {n — 1) dimensions, and irrel- 
evant for d > 2/{n — 1) Hence we identify the up- 
per critical dimension, below which fluctuations in fact 
dominate the asymptotic behavior, as dJn) = 2/{n — 1) 
for n-th order annihilation processes pl| , p^ . Thus, for 
n > 3 fluctuations are not too important in any physical 
dimension d > 1. ^ 

Furthermore we notice that the coupling D is irrele- 
vant, i.e., compared to the other parameters in the theory 
its influence should become negligible in the asymptotic 
long-time, long- wavelength limit. Evidently, pB(x,t) ei- 
ther vanishes (inactive phase) or approaches a constant 
(active phase) as i ^ oo. In the former case, normal 
B diffusion, if present {D' > 0), will dominate; in the 
latter situation, the combined quantity D p"^ will effec- 
tively act as an ordinary diffusion constant, numerically 
renormalizing D' . In any case, we see that the ordinary 
A diffusion process is not qualitatively affected by the 
induced hopping through attractive coupling to the B 
density, and the associated noise cross-correlations. Also 
when D' = 0, as in our original model, and in a system 
with an initially finite number of B particles, asymptot- 
ically the A particles either remain localized or display 
standard diffusion. In this respect, in numerical simu- 
lations the induced anomalous diffusion in which we are 
interested here would appear as a crossover feature in 
the long-time kinetics and correspond to corrections to 
scaling to the leading asymptotic time dependence. In an 
infinite system, however, with initially finite B density, 
the anomalous diffusion regime will persist indefinitely. 

A corollary of these observations is that the rate D 
does not acquire any non-trivial frequency or time de- 
pendence in the infrared. In the field theory language, 
we note that neither diffusive propagator for the A or 
B species can be renormalize d by t he {abba) four-point 
vertex in the unshifted action ( ^.13 ) , or equivalently, the 



three- and four-point vertices in the shifted action ([3.l|). 
Consequently, the renormalization for the vertex func- 
tions T-gj^ or r^g^^ and Taba, respectively, can be deter- 
mined to all orders in the perturbation expansion (with 
respect to D) by means of a Bethe-Salpeter equation, or 
equivalently, a geometric series of loops containing just 
the A and B propagator. This leads to the renormalized 
wavevector- and frequency-dependent coupling 



DR{q,uj)^D 



D 



d'^p [((f — p)'^ 



P 



{2n) 



(3.6) 



where q and uj denote the momentum and frequency 
transfer between the A and B particles. We may now 
set D' — Q again, and investigate the long-wavelength 
limit (f — > 0, 



dq 



^DR{q,uj) 



= D 



9=0 



D /UJ 
D \D 



d'^k fc2 
{2tj:Y -i + P 



(3.7) 



where k"^ = Dp^/uj. Thus, as ^ 0, the fluctuation 
corrections vanish (provided the integral is regularized 
in the ultraviolet with an appropriate cutoff), and the 
renormalized coefficient Dr in Eq. ( |3.7| ) approaches the 
original "bare" constant D. 

This is to be contrasted with the infrared-singular be- 
havior of e.g. the B pair annihilation rate, for which an 
analogous procedure yields ||l8|] 



Aj?(0,w 



= A 



1 + 4 



d'^p 



1 



H -1 



D J {2ttY 



D 



-iuj/D + q^ /A+p'^ 

-1 



1 



d'^k 



{2ttY -i + 



(3.8) 



For d > dc{2) = 2, again Afl(0,0) = A is just the origi- 
nal rate constant, resulting in the mean- field power law 
PB{t) oc t^^. However, for d < dc{2) = 2, Xii{0,uj) cx 
(^1-^/2 vanishes for low frequencies. Inserting the corre- 
sponding effective time-dependent rate XR{t) cx t-i+^^Z^ 
into Eq. (3.3) leads to the correct slower algebraic decay 

PB{t) (X t^\ 

In summary, the B process itself is, per definition of 
our model, not influenced by the A dynamics. In the 
renormalization group treatment, this is reflected by the 
fact that the coupling D is irrelevant, and thus does not 
affect the long-time behavior. Yet the induced hopping 
rate D is of course crucial for the A species kinetics, and 
must be kept even in the mean-field approximation. We 
may thus solve for the B kinetics first, and then explore 
its influence on the induced A diffusion. Henceforth, we 
shall again set D' = 0, as otherwise simple ordinary A 
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diffusion would ensue, with D then irrelevant also for the 
A kinetics, and the entire coupling of the A and B pro- 
cesses would disappear asymptotically. In the following, 
we shall study the A kinetics assuming a spatially homo- 
geneous, but time-dependent distribution of B particles, 
which leads us to a mean-field description. 



IV. MEAN-FIELD EVOLUTION EQUATIONS 

A. Annihilation kinetics 

Let us assume we can ignore spatial fluctuations for 
the B species entirely, and ignore the reaction noise. 
For the n-th order annihilation processes, we saw that 
this as at least a qualitatively correct description for 
d > dc{n) = 2/(n — 1), i.e., for d > 2 in the case of pair 
annihilations, d > 1 for the third-order process ?>B — > 0, 
and in any physical dimension for n > 3. The evolution 
equation can either be obtaine d di rectly from the non- 
equilibrium operator in Eq. ( |2.3| ) and the equation of 
motion (2.9), or from solving for the stationarity condi- 
tion 5AI5 h = for the action A = Ab+ Ar, Eqs. ( |2.10| ) 
and (|2.1l| ), setting L» = 0. (Notice that 5A/6h = is 
always solved by & = 1.) Either procedure results in the 
obvious mean-field rate equation 



dt PB{t) = -n\n psitT > 
which is readily integrated for n > 1, 



PB{t) 



Pb(0) 



Ps(0)i 



(l + i/r)'/^"-'^ ' n(n-l)A„' 



(4.1) 



(4.2) 



i.e., for t ^ r the B density decays algebraically cx 
in this approximation, while of course for n — 1 



PB{t)= PB{0)e-^'' . 



(4.3) 



In the same manner, we may obtain the evol ution equa- 
tion for the A species, or just consider Eq. ( |3.2| ) for D' — 
and vanishing noise. In the spirit of mean-field theory, 
we assume a homogeneous B density, and obtain 



dtPA{x,t) = DpB{t)W^pA{x,t) 



(4.4) 



Again, this equation can be solved exactly, considering 
a delta-like density distribution for the A species at the 
initial time t — 0. As in this mean- field approach the 
PB{t) is assumed to be spatially uniform, the A species 
will be Gaussian distributed in space, just like in ordinary 
diffusion. 



Pa{x, t) 



d/2 



exp 



2^{xl{t))J 2(fi(t)) 



(4.5) 



However, the B decay (or lattice depletion) will be re- 
flected in the anomalous time-dependence of the width 



(mean-square displacement). A straightforward brief cal- 
culation yields 



{x\{t))^2D / pB{t')dt' 



(4.6) 



For n = 1, i.e., the simple exponential decay (4.3), the 
result is 



(4(0) = ^(1 



(4.7) 



Initially (Aii ^ 1) one finds normal diffusion with ef- 
fective diffusion constant D — D pb{Q), but at long 
times the mean-square displacement approaches a con- 
stant, and the A particles remain localized in a region of 
volume oc {x\{t oo))'*/^ ^ {2D pb{0)/ Xi f^^ . Given 
that this simple process is characterized by short-range 
correlations in space and time only, we do not expect any 
considerable modification through fluctuation effects. 
In the pair annihilation case, n = 2, one finds 



{x\{t)) = 2DpB{0) In 1 + 



t 

T 

while the mean-field result for n > 2 reads 



~ n — 1 

{x\{t))^2DpB{0)^^T 



(4.8) 



(4.9) 



In the asymptotic regime t ^ t, this implies anomalous 
diffusion according to 



(4(t)) cxt2/(2+e) 



(4.10) 



with a positive exponent O ~ 2/(n — 2) indicating sub- 
diffusive behavior. In the limit n — > oo we have — > 0, 
and conventional diffusion is recovered. The reason is of 
course that for large n the depleting reactions become 
very unlikely, as n particles are required to meet at the 
same lattice site. Thus, low-order B species reactions 
are much more effective in slowing down the A diffusion. 
The time scale for the crossover to the pure algebraic 
decay of the B particle density and subsequently for the 
anomalous A diffusion is given by r oc /9b(0)^~"/A„. The 
crossover to the asymptotic slow dynamics is fast for large 
initial densities and reaction rates. 

The above analysis should be qualitatively correct for 
n > 3, as the corresponding critical dimension dc(n) < 1. 
For n = 2, i.e., B pair annihilation processes in d < 2 
dimensions, we know that at long times anficorrelations 
develop j2^jl^,|l^: Initially close-by particles disappear 
quickly, and only widely separated ones survive. This 
effective "repulsion" should result in a roughly uniform 
spatial B distribution even for a clustered initial config- 
uration. Given that the coupling coefficient D itself does 
not renormalize, we therefore expect that our decoupling 
assumption leading to Eq. (4.4) should represent a fair 
approximation, provided the correct time dependence of 
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the B density is inserted. For d < 2, the asymptotic 
result is 



PB{t) OC t-"^^ 



(4.11) 



see Ref. [|l8|| and also Sec. Ill following Eq. (3.8), whence 
{x\{t)) = 2Dt^-'^/^ (4.12) 

with an appropriate effective rate D (x D/{1 — d/2). In 
low dimensions, this algebraic subdiffusive behavior with 
Q = 2d/{2 — d) replaces the logarithmic law (4.?). At 
the critical dimensions dA2) = 2, one finds the typical 
logarithmic corrections |18|] 



implying 



PB{t) (X t ^ Int 



(4.13) 



(4.14) 



which also descri bes slower kinetics than given by the 
mean-field result (4.6). For the case of n = 3 at its criti- 
cal dimension dc{3) = 1, 



PB 



(t) cx {t-^ lni)i/2 



(4.15) 



and one would therefore expect the leading time depen- 
dence 



{x\{t)) oiDit lnt)i/2 , 



(4.16) 



i.e., essentially a square-root power law with logarithmic 
corrections. 



B. BARW kinetics 

We now extend the B dynamics and include branch- 
ing processes of the for m B — > (m + 1) B with rate am, 



described by Eqs. (2.4) or (2.12). The mean-field rate 
equation ([1.1|), with rt > 2, is then replaced by 



dtPB{t) = -nXnPBity +ma„iPB{t) 



which has two stationary solutions pb — ^ (inactive 
phase) and 



Pb = 



l/(n-l) 



(4.18) 



(active phase). For any am > 0, the latter turns out to 
be stable, i.e., BARW are always in the active phase in 
the mean -field approximation. The explicit solution of 
Eq. (4.18) furthermore shows that the asymptotic den- 
sity is exponentially approached. 



Psit) 



oo 

Pb 



[l + Cf 



-(n-l)rnCT„tll/("-l) ' 



(4.19) 



where C = [/og'/pi3(0)]""^ - 1. Again, Eq. (jij) is solved 
by the Gau ssian distribution (4^) with mean-square dis- 
placement ( f4.6D . The ensuing integral is readily calcu- 
lated for some special cases, e.g., for n = 2 



{xl{t))^2Dp^ 
whereas for n — S 

bpi 



t+ ln(l + Ce-""""*) 



= 



In 



Vl + Ce- 



1 



\/l + Ce-2""^"*- 1 



(4.20) 



(4.21) 



In general, asymptotically normal diffusion with effective 
diffusion coefficient D p"^ is recovered in the active state. 



{x\[t))^2Dp^t 



(4.22) 



The properties of the active phase with an asymptoti- 
cally homogeneous B density are n ot much influenced by 
fluctuations, and hence Eq. ( 4.22 ) should aptly describe 
the ensuing A kinetics even beyond mean-field theory. 

For the possible existence of an inactive phase, and the 
characterization of the ensuing critical behavior, fluctua- 
tion effects arc however of utmost importance for n — 2, 
and it turns out that the cases of odd and even offspring 
number m need to be distinguished. For odd m, aside 
from all lower-order branchings, first-order decay pro- 
cesses B % are generated, and become sufficiently 
efficient to shift the critical point to CTc > for d < 2 
dimensions. The emerging transition at CTc > can be 
shown to be in the generic directed-percolation (DP) uni- 
versality class The inactive phase is then governed 
by exponential B density decay, whereu pon the A species 
will become localized according to Eq. (|4.7| ) . At the crit- 
ical point itself, the B species density decays according 
to a power law pb{S, t) ~ with a — (ij zv^ given by 
DP critical exponents in = 1 and c? = 2, respectively. 



(4.17) This would suggest {x\(t)) oc yet the B density 



is far from uniform at the critical point, and is instead 
characterized by the appearance of fractal density clus- 
ters. While we would still expect subdiffusive behavior 
for the A species with 8 > 0, this exponent will likely be 
influenced by the power-law correlations in the critical B 
density. In the case of even m, for which the B particle 
number parity is locally conserved under the reactions, 
a non-trivial transition with ctc > is possible only for 
d < « 4/3 dimensions. The inactive phase is then 
given by the pure pair annihilation theory, and conse- 
quently Eq. (4.12) should provide a fair description for 
the ensuing anomalous A diffusion. The critical behavior 
is governed by a different parity-conserving universality 
class, with a < 1/2. In this instance, we again expect 
the above mean-field description to be rather inaccurate. 
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V. CONCLUSIONS 



We have studied a novel mechanism to induce anoma- 
lous diffusion. Whenever an active particle of the A 
species performs a random walk on a lattice, it may 
visit a certain lattice site only provided this site is al- 
ready occupied by at least one B particle. The random 
walk is prolongated when the B particles react with each 
other in such a manner that the B species density is 
decreasing. If that decay is exponential (first-order re- 
action), then after a short time interval (given by the 
inverse decay rate) the B species has disappeared and 
a further visit of an A particle at that site is impossi- 
ble. As a consequence the A species, after some initial 
mixing, remains localized. When the B species undergo 
reactions of higher order, requiring at least two B parti- 
cles to meet at a lattice site, an algebraic decay ensues 
that allows hopping processes for the A species to occur 
for a much longer period. However, the random walk 
process is slowed down considerably as the B density 
diminishes, resulting in a much shorter mean-square dis- 
placement of A particles as compared with conventional 
diffusion. The emerging anomalous diffusion is governed 
by power laws or logarithmic behavior that can (approx- 
imately) be related to the asymptotic time behavior of 
the reacting B particle density. In this instance one may 
view this process as resembling diffusion on a dynamical 
fractal. Only when at long times the B density remains 
finite and nearly homogeneous, conventional A diffusion 
is recovered. This situation corresponds to diffusion with 
dynamical disorder, where in the long-time limit the B 
particles, with largely decayed fluctuations, merely re- 
semble a quasi-static inhomogeneous background for the 
A kinetics. The consistent mathematical treatment of 
diffusion on a static fractal, as well as induced diffusion 
processes on critical (isotropic or directed) percolation 
clusters or near BARW critical points remains an open 
problem that requires more sophisticated analysis beyond 
the largely mean-field approach presented here. 
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